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Abstract 

C^ . We obtain the conductance of a system of electrons connected to leads, within time-dependent density- 

functional theory, using a direct relation between the conductance and the density response function. Cor- 
C I rections to the non-interacting conductance appear as a consequence of the functional form of the exchange- 

^ ■ correlation kernel at small frequencies and wavevectors. The simple adiabatic local-density approximation 

and non-local density-terms in the kernel both give rise to significant corrections in general. In the homoge- 
> 
C^^ ' neous electron gas, the former correction remains significant, and leads to a failure of linear-response theory 

^ ■ for densities below a critical value. 

o 
\o 
o 
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Time-dependent density-functional theory extends the domain of ab-initio calculations to sys- 
tems carrying a current, but relies on the accuracy of the exchange-correlation energy functional of 
the electron density (at present and past times), which in practice has to be approximated. Impres- 
sive success has been achieved within the non-equilibrium Green's function formulation using the 
simple ground-state density-functional exchange-correlation potential in a self-consistent formula- 
tiofti^ (gDFT). However, limitations of the latter approximation were recently identifie d^i'^1^1^1^1^ . 
For instance, gDFT's omission of the derivative-discontinuity in the exchange-correlation energy 
functional was found responsible for serious errors in transport calculation through localized res- 
onant levels^. Improvements through an unrestricted gDFT formulation have been argued to 
describe properly some aspects of the Coulomb blockade in quantum junctions^. At this level of 
the theory the exchange-correlation potential of the equilibrium system, v^c, is responsible for the 
electron interaction effects. 

In a further theoretical development, Na Sai et ah— identified a dynamical correction to the re- 
sistance of a quantum junction stemming from the contribution of the exchange-correlation elec- 
tric field to the overall drop in the total potential, as reflected in the exchange-correlation kernel 
fxc- They estimated the correction within time-dependent current-density functional theory^ (TD- 
CDFT) and showed that it has its origin in the non-local density-dependence of the functional. 
The very applicability of time-dependent density-functional theory (TDDFT) to the problem of 
quantum transport in the long-time limit has been discussed in depth by G. Stefanucci and C.-O. 
Almbladh^ and by M. Di Ventra and T. Todorov^. 

Several authors have proposed alternative treatments that avoid the complexities of the 
exchange-correlation kernels of TD(C)DFT, either by using the usual gDFT approach in com- 
bination with a model self-energy within the central region ^*^, or by treating the central region 
with the configuration integration method^ ^ while approximating the non-equilibrium distribution 
of the electrons. 

In this paper we address the class of corrections to the linear-response conductance that arise 
from the exchange-correlation kernel and potential in TDDFT. A kernel must give a satisfactory 
description of the linear-response regime if it is to be of use in more general quantum conduc- 
tance calculations. Apart from the ultra-non-local contribution found by Sai et al?, we identify a 
new correction to the conductance that already appears within the adiabatic local-density approx- 
imation (ALDA), and gives a significant increase in the conductance even for the homogeneous 
electron gas. By using an expression for the conductance in real space, we can explore both the 



ALDA correction and the correction of Sai et al. for homogeneous and inhomogeneous systems, 
in relation to the gDFT conductance. 

The central concept in our approacbi^ is the identification of the conductance G as the strength 
of the Drude singularity of the conductivity tensor in reciprocal spaced, 

\ima,,{q,q',Y.^^0-(o)^lJlG5{q)5{q'), (1) 

where we consider a geometry where z is the direction of the current flow, q a reciprocal vector 
in that direction, and for simplicity we assume that the system is translationally invariant along 
the x,y directions and Kj^ is a wavevector reciprocal to R^ = {^,y) , i-C- we consider an ideal 
interface. The conductance G is then a conductance per unit area of the interface. The limit 
(0^0 must be taken from the upper half of the complex frequency plane ft) = 9l{ft)} + ?'a, a > 
as the last step in the calculation. This order of limits - first an infinitely long system (i.e. q 
continuous), and only then ft) ^ - is essential, for otherwise different and even divergent results 
are obtained, arising from the incorrect "piling-up" of charge at the ends of the system^^. The 
unambiguous evaluation of a finite conductance for an infinite dissipationless system is facilitated 
by the adiabatic switching-on of an external electric field, characterized by a finite drop in potential 
/^yext^ that can convenientl)*^^ be represented as the field 

Ayext j^ 

Ef{z.t) = -fi^^(^e-'«^ (2) 

for the time interval t G (— °°, 0). c is a constant with units of velocity and a > controls the speed 
of switching on the field. 

Through Eq. (HJl we can relate the conductance to the non-local conductivity for small ft). In 
turn, the conductivity is simply related to the irreducible polarizatio n ^ ^i ^ 3 26 

O^{q,q',K^=Q-0)) = '-^t'\q.q'-^Gi)- (3) 

qq' 

The polarization is conveniently calculated via the density response function calculation within 
TDDFT (although other possibilities exist, such as TDCDFT or many-body perturbation theory). 
The polarization function satisfies the Dyson equation 

t''\q,q';G))=X^{q.q'\Gi) + 
Jdq"dq"'x\q,q";CO)fUq",q"';CO)x"'''iq"',q';CO), (4) 



where X^{<1j1'''>(0) is the non-interacting Kohn-Sham density response function. The above two 
equations give the conductance of an interacting electronic system: for a given kernel f^c one 
needs to invert the Dyson equation Q and substitute the result into Eq. dSj. 

However, we can gain insight by multiplying Eq. Q by ico/{qq') and taking the limit CO ^ 0. 
Clearly, the resulting left-hand- side is singular in q, with the strength being directly the conduc- 
tance of the interacting system, G. The strength of the first term on the right-hand- side of dH 
multiplied by the same factor gives the conductance of the non-interacting Kohn-Sham system, 
G^. The difference between these two conductances, i.e. the correction due to the exchange- 
correlation kernel, is then nonzero only if the last term in Eq. ^ also leads to a singular form. 
This observation can be used to deduce the aspects of the kernel that do influence the conductance, 
since the general character of x for small O) is well known. 

The most obvious choice, making use of the character of x^ '' ~ 'JQ'^ is fxc {q^q'^O}) ~ 
-^A(q,q'; CO) where A{q,q'; co) ^ A = A{q = 0, <?' = 0) 7^ for ft) ^ 0. The resulting con- 

qq 

ductance then takes the form 

From the above expression we see that A represents some part of the dynamical resistance per 
unit area and, in fact, it is equivalent to the correction identified by Na Sai et al.—. It can be 
shown ^^ that a purely longitudinal exchange-correlation electric field used in their treatment^ii^ 
within TDCDFT is equivalent to a contribution to the TDDFT kernel of the asymptotic form for 
small q, q' 

qq' J 3 V '^UJ / W 

where n{z) is the number of electrons per unit cross-sectional area and 77 is the dynamical viscosity 
of a homogeneous electron gas^. We should note that for homogeneous systems A = since 
dyn{z) = 0. This is important since the functional form fx" given above and the limiting process 
would not lead to a finite result for a homogeneous system. 

However, the form above does not exhaust all the possibilities. Most surprisingly, the con- 
ductance is also affected by a local-density term f^c {q,q';(o) = f^c {q — q'\Oi) -^ B{(o)d{q — 
q'),B{o}) — > 5 7^ for ft) ^ 0, which, Fourier-transforming f^c to real space, naturally appears 



within the adiabatic LDA^^ in an extended system, 

fi\v,v')=fi^''^[n\v)]5{T-v') (7) 

f^^^^%^]=d\n'^e,c{n^))/dn\ (8) 

(Here e^c {n^) is the exchange-correlation energy per particle of a homogeneous electron gas (HEG) 
of density rP). To exhibit this case let us consider the ID HEG where the algebra becomes partic- 
ularly simple, and we will assume that the kernel B\o = /^^[n'^] exists even for this non-Fermi- 
liquid system. (For 2D and 3D cases, where this fact is generally accepted, the analogous steps 
lead to technically more demanding formulas and only some of the results can be obtained analyt- 
ically. However, we will demonstrate numerically that the qualitative general picture is identical 
to that described below for the ID situation.) 

The non-interacting response of the ID gas, Xij)^ for q «kf, has the simple form 

When fxc is combined with this in the Dyson equation ©, a renormalization of the conductivity 



results, 

Oij){q,q'-C0)^^- ^B,^^',"^ ^ -x5{q-q'), (10) 

which, using Eq. dU), gives the conductance 



2k f —ico 



GiD= .^°° (11) 



Ttkf 

Typically (i.e. in 2D and 3D) the adiabatic kernel (EJ) is negative and decreasing function of the 
density and eventually as a certain critical density (characterized by kf^c given by ^5id (^f.c) ~ 
— 1) is approached, the conductance shows a singular increase. This transition corresponds to the 
known instability of the HEG against arbitrarily small fluctuations in the total potentiaU^, signified 
by the appearance of a pole in the upper half of the complex plane of the irreducible density 
response x"^"^- For densities beyond the transition, linear-response theory is therefore inapplicable. 
The situation in 3D gas is qualitatively similar. For q << qp the non-interacting response has 
the form (ft) = ia, a > 0) 



0.03 



e) 




FIG. 1: (Color-online) Dependence of the conductance per area of the Kohn-Sham gas in 3D, G^, and the 
real and imaginary parts of the full conductance within the ALDA, {G^^'), on the density parameter, r^. The 
ALDA correction is significant within the physically relevant region 2 < r^ < r^ c- 



from which, by means of Eqs. ([T]) and Q, we obtain the conductance per unit area of the Kohn- 

n k^ 
Sham gas G = j^. The critical density at which the instability occurs can be found from the 

appearance of poles of the the irreducible response function in the upper half of the complex plane, 

2 

which leads to the criterion B{kp^c) < ~j~- The form of the correction to the conductance in 3D 
cannot be obtained analytically, so we show in Figure [T] our numerical results, which qualitatively 
resemble the ID case. We see that the ALDA correction leads to a systematic increase in the 
conductance. We stress that this correction is a direct consequence of the proper order of limits 
performed in Eq. |4l From this it also follows that this correction would not be present within the 
static gDFT calculations based on the NEGF formulation but would appear within a direct time- 
evolution approacbi^ii^ if the ALDA is employed. (A minor difference in 3D is that beyond the 
critical rs^c the formally defined conductance retains a nonzero real part, whereas in ID it is pure 
imaginary.) 

We briefly discuss a third form of the exchange-correlation kernel, f^c {q,q'',CO) ~ 

2 

^C{q,(o)5{q — q'), C{q,0)) ^ C 7^ for ft),^ ^ 0. It can be shown that this form of fre- 
quency dependence in the kernel also changes the final conductance. This kernel is similar to the 
one used by Botti et al. ^^ and Reining et al?^ for bulk insulators and semiconductors characterized 
by a finite gap in the electronic density of states at the Fermi energy. In their estimate, which lead 
to significant improvement in the electron energy-loss spectra and optical absorption spectra, the 
coefficient C ~ ft)^ where (Og is the effective energy gap. This suggests that fxJ = will be zero, 
or small, in a metallic system. 



In order to explore the relevance of these conductance corrections for inhomogeneous systems 
we consider the metal- vacuum-metal case of two jellium surfaces separated by a distance d. In this 
case the parameter B has been evaluated using Eq. ^ and the Perdew-Zunger^^ parametrization 
of the quantum Monte Carlo correlation energy of the HEG^^. The parameter A is obtained from 
Eq. ®, and its contribution to the conductance follows from Eq. Q- 

In our calculations we employ two jellium slabs of thickness L and density given by r^. The 
calculation of X^{z:z!\icc) is performed at the self-consistent LDA level ^^. Subsequently we 
invert the Dyson equation, Eq. ©, in real space and thereby calculate the irreducible response 
X"^{z, z!\ ice) . For the final step we employ an integral form of Eq. dH) 

rdqdq'x'''''{q,q';icc) 
G=- hm a — , (13) 

a^0+ J 2n qq 

where we have used the identity ^. Direct implementation of this expression in reciprocal space is 
numerically rather cumbersome, since x"^ is, in practice, calculated for a finite (periodic) supercell 
for which the region with q,q' « I h poorly described. This problem can be circumvented by 
re-expressing G in a real-space form, 

G=-\ima dz dz!x'"''{z,z!\ia), (14) 

which is obtained from (IT3t by utilizing the fact that the 1 /qq' singularity is only apparent since 
X{q-,4) ~ w' for small q,cl . Extrapolation to zero frequency (from w ~ 0.01 Ha << Ep^ is 
done in parallel with extrapolation of the thickness of the slabs L to infinity. Further details of 
the calculation will be published elsewhere^. The quality of the numerical procedure can be 
judged from the correct exponential decrease in the Kohn-Sham conductance over several orders 
of magnitude shown in Fig. El 

The resulting dependence of the conductance on the vacuum width is shown in Fig. |2l for a 
representative jellium density r, = 3 a.u. (Au). For small vacuum widths J, the correction due 
to the ALDA kernel f^c clearly dominates; for larger widths (J ~ 6 — 8 a.u.) /i" and fxc shift 
the Kohn-Sham conductance in opposite directions and to some extent cancel each other. We 
should note that unlike the use of a global value for the viscosity in the original work by Sai^, we 
have used the local viscosity r\ ^f{z)\ in Eq. Q, determined by the local density, which leads to 
suppression of the non-local corrections for d > 2a.u.— 

In conclusion, we have presented a unified formalism, based on the singular character of re- 
sponse functions, to address the conductance of a general system of interacting electrons. We have 
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FIG. 2: (Color-online) The Kohn-Sham conductance G° (upper panel) and the corrected conductances G 
(shown relative to G°) (lower panel), as a function of the vacuum width d. The corrections due to the 
non-local (a) and ALDA (b) kernels have opposing signs. 

explicitly identified three different contributions to the dynamical resistance: (a) the non-local con- 
tribution^ parametrized by the dynamical viscosity of the homogeneous electron gas, effective only 
for inhomogeneous systems, (b) a local contribution parametrized by the adiabatic LDA exchange- 
correlation kernel, and (c) an ultra- non-local contribution that, from presently available estimates, 
is not important for metallic systems. These three forms are unlikely to exhaust all the possibilities, 
and our theoretical framework remains applicable for further analysis based on many-body pertur- 
bation theory to obtain additional relevant contributions to the exchange-correlation kernel or local 
field factor. In the homogeneous electron gas, we have found that the conductance diverges as the 
critical density rs^c is approached, beyond which the gas is unstable. For practical calculations, 
our formalism can readily benefit from available computer codes for ab-initio calculations in real 
materials that can generate density-density response functions. We have shown the importance of 
the non-vanishing corrections for the conductance of a model inhomogeneous system. 
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